Amplification of realistic Schrodinger cat-like states by homodyne heralding 
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We present a scheme for the amphfication of Schrodinger cats that collapses two smaller states 
onto their constructive interference via a homodyne projection. We analyze the performance of the 
amplification in terms of fidelity and success rate when the input consists of either exact coherent 
state superpositions or of photon-subtracted squeezed vacua. The impact of imprecise homodyne 
detection and of impure squeezing is quantified. We also assess the scalability of iterated amplifica- 
tions. 



Coherent state superpositions, or optical Schrodinger 
cat states, are widely recognized as promising resources 
in quantum information pJ-l5J, quantum metrology [6^- 
[8j, and fundamental tests p1-[T2l. In the near-orthogonal 
basis of coherent states (7|— 7) = e~^^ , two particular 
instances for these states are 

l^±(7)) = , \, „ (|7) ± 1-7)) , (1) 

where the sign (±) of the superposition refers to the 
even and odd cat state, respectively. These states ex- 
hibit quasi-probability distributions in phase space which 
are distinctly non-classical. This makes them all the 
more challenging to generate deterministically as that 
would require strong Kerr- type non-linearities [TSHlS]. 
One has then to resort to heralding techniques which, 
though probabilistic, need only linear optics and projec- 
tive measurements [16^ . These state-engineering schemes 
are nonetheless approximative and present a limita- 
tion in the fidelity they produce with ideal cat states. 
Photon-subtraction of squeezed vacuum, for example, is 
a well-established method to generate approximations of 
small amplitude cat states, colloquially referred to as 
Schrodinger kittens p!7H2T] . Even in the best experimen- 
tal conditions, the fidelity between the photon-subtracted 
squeezed vacuum (PSSV) and an actual odd cat state 
|/^_(7)) degrades markedly for 7 > 1.2 [22 . Yet, for 
these states to be reliable resources in quantum compu- 
tation, their fidelity with cat states at least as large as 
7 = 1.2 need to be maintained at near-unit fidelity [3j[23]. 

Single-photon subtraction is only one example of sev- 
eral measurement-induced schemes which have been pro- 
posed to generate kitten states [20l I24H29] . However, 
none of these schemes can produce arbitrarily large cats 
in a single run. Ways to get around this issue have been 
devised using the recursive amplification of small, ap- 
proximate cats [30, 31 . For example, it was suggested in 
[32] to interfere a supply of delocalized single photons fol- 
lowed by homodyne heralding to generate large entangled 
cat states. These proposals have in common that they 
rely on the coherent mixing of two small cats, where- 
upon a projective measurement collapses one of the two 



outputs onto a constructive interference of the inputs — 
hence the amplification. Here, we shall pursue the same 
idea but make use solely of homodyne heralding for its 
relative simplicity and high quantum efficiency. We also 
demonstrate that the acceptance window of homodyne 
heralding can be widened to increase the success rate of 
the amplification while at the same time maintaining a 
satisfactory fidelity at the output. 

The outline of this article is as follows. In Sec. Uwe 
review the generation of odd kitten states from squeezed 
vacuum. The output is compared to the ideal odd cat 
state and the effects of impure squeezing are illustrated. 
In Sec. [lljwe present the amplification scheme in the case 
of ideal input cats and model the effect of a wide homo- 



dyning window. Sec. [III|then considers the amplification 
of the more realistic PSSV. The impact of both impure 
squeezing and wide post-selection is illustrated. In [TVl 
we return to the idealized case of ideal homodyning and 
pure squeezing to consider how our scheme scales with 
iterated amplifications. 



I. APPROXIMATION OF SMALL ODD CATS 

In this section, we shall briefly review the genera- 
tion of odd kitten states from the photon-subtraction of 
squeezed vacuum and analyze its performance in the face 
of impure squeezing. The basic setup is depicted in Fig. 
[1] The original proposal of Dakna et al. [17] required 
that the photon subtraction be performed by photon- 
number resolving detectors. However, as is done in most 
practical schemes [18l [TQ] [JT , we shall assume that the 
post-selection is a binary detection of either zero pho- 
tons or at least one photon, as would be allowed for by 
an avalanche photodiode (APD). The modeling of such 
an "on-off" post-selection operation is presented in Sec. 
lA3l 

Quantum inefficiencies and dark counts are two nui- 
sances inherent to photodetection which should be reck- 
oned with. Whereas the former merely affects the suc- 
cess probability of the scheme, the latter contaminates 
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the output with a squeezed vacuum component which 
weakens the nonclassicahty of the output. An equahy 
detrimental effect is the impurity of the squeezed vac- 
uum. We shah not treat quantum inefficiencies and dark 
counts here as they have aheady been covered in [16 in 
the context of cat amphfication. We wih however look 
at the fidelity response to the amplitude and impurity of 
squeezing. 



APD^ 

/post-selection 



|0> 



FIG. 1: (Color online) Setup for the generation of approxi- 
mate small odd cat states. A squeezed vacuum, represented 
here by the pumping of a nonlinear medium, is partially re- 
flected onto an "on-off " photon detector such as an avalanche 
photodiode (APD). Upon the reflection of at least one photon 
from the squeezed vacuum, a state that likens the squeezed 
single photon is prepared in the limit of unit transmission 
1. 
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A. Fidelity of the output for pure squeezing 

Let's assume from now on that the squeezed vacuum 
is pure and anti-squeezed in the x-quadrature. I.e., if we 
denote the squeezing relative to the shot noise variance 
by ^ (in dB), then — — > 0, where the subscript 
labels the measured quadrature [38^. The fidelity of the 
PSSV state with an ideal cat state \K-(a)) then depends 
on the squeezing ^ and the transmission T of the subtrac- 
tion beam splitter [39 . The effect of these parameters is 
illustrated in Fig. [2] where we can see that the fidelity 
is optimized for low squeezing and near-unit beam split- 
ter transmission T ^ 1. For any given input squeezing, 
there corresponds a finite amplitude a of the target cat 
with which the output has an maximized fidelity. For 
example, a squeezing of, say, 3 dB is optimal for produc- 
ing an approximation to a cat state of size a = 1.0. A 
complementary investigation of PSSV states that looked 
at nonclassicahty as the main figure of merit (as opposed 
to fidelity) is given by Kim et al in [33] . 

In order to assess the performance of the PSSV gen- 
eration, it is most instructive to restrict ourselves to the 
maximal fidelity achievable with any target cat of size 
a. This fidelity optimum is traced by the ridge of Fig. 
[2] (top) and its a dependence is reproduced in Fig. |2] 
(bottom) along with the corresponding success rate. For 
example, if one wants to produce an approximate odd cat 
state of size a = 1.5, the required squeezing should be 
^ ~ 5.2 dB for a fidelity of at most 95.4% and a success 
probability of 0.4%. 

As far as the success probability is concerned, it can 
be increased at the expense of fidelity by increasing the 



FIG. 2: (Color online) Top: Contour plot of the fidelity of 
the PSSV state with \Hi-{a)) given an input squeezing be- 
tween and 15 dB. Bottom: Maximum achievable fidelity 
of the PSSV with \Hi-{a)) (black) and corresponding success 
probability as a function of a (gray). In both cases, the trans- 
mission of the subtraction beam splitter is 99%. 

incidence of photons on the APD via stronger squeezing 
or weaker beam splitter transmission. 

B. Impact of impure squeezing 

In any real world experiment, noise will inevitably 
infiate the variance of the squeezed vacuum in either 
quadrature. This impurity has been explained as stem- 
ming from losses or from multimode parametric down- 
conversion [34 whereby the photons triggering the post- 
selection belong to a different spatial or frequency mode 
than the heralded state. Impure squeezing may origi- 
nate in the down- converter itself or, more generally, at 
any point in the setup where vacuum contamination or 
modal mismatch could take place, including at the detec- 
tors (e.g., via quantum inefficiencies). Regardless of its 
root causes, we shall wrap these impurities into a single 
parameter e relating the squeezing ^ in dB of the x- and 
^-quadratures, 

= -e^p, (2) 

whereby pure squeezing corresponds to e = 1. As derived 
in §A 1[ the purity of the squeezed vacuum is given by 

V = IQ-^^^-'^^p. (3) 
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Note that the Heisenberg uncertainty relation imposes 
that e > 1. 

The impact of impurity on fidehty is plotted in Fig. |3] 
where we set the squeezing at —3.0 dB and adjust the 
anti-squeezing according to four different settings of pu- 
rity. The immediate observation is that a decrease of 
10% in purity, from 100% to 90% — corresponding to an 
increase of anti-squeezing to +3.9 dB from +3.0 dB, leads 
to a drop of the maximum fidelity of nearly 32%. 
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FIG. 3: (Color online) Fidelity between an ideal cat state 
and the kitten states obtained from squeezed vac- 
uum = — 3 dB at four squeezing purities V (marked in 
decreasing darkness of gray): 100%, 90%, 80% , and 70%. 
The Wigner functions of the PSSVs corresponding to these 
four purities are shown at the top. As a reference of what 
the "ideal" output ought to look like, the Wigner function of 
the state \k,-{1.06)), which has the highest fidelity with the 
PSSV obtained from pure squeezing, is shown at the bottom. 
The transmission of the subtraction beams splitter is set at 
99%. 



II. AMPLIFICATION OF IDEAL ODD CATS 

In Fig. [4] we present our amplification setup: Two iden- 
tical cats are mixed on a balanced beam splitter where- 
upon one of the ensuing modes heralds the amplified out- 
put based on the measurement of x-quadratures around 
X = 0. The scenario where the inputs are idealized cats of 
opposite parity has been outlined by Takeoka and Sasaki 
in [31,. We will however look at the more practical case 
where the input have identical parity. 

Let's briefiy run through the evolution of the state in 
setup. We can readily see that the state emerging from 
the first balanced beam splitter contains an even cat of 
amplitude a/2 times larger: 



\i^±{a)) (g) \f^±{a)) 

\tz+{V2a)) |0) ± |0) \i^+{V2a)). (4) 



If, by measuring one of the modes, we can post-select 
|0) from |/^+(>/2a)), then the other mode will collapse 
onto the desired amplified state \tz-^{y/2a)). The accu- 
racy of this discrimination is of course limited by the 
overlap of vacuum with the cat {0\K-^{^/2a))^ which is 
however negligible for a ^ 0. Based on the wave func- 
tion profiles of the vacuum and the cat, we can see that 
homodyne measurement of the x = quadrature is in- 
deed a good way to discriminate the two states as it is 
where their overlap is minimized (cf. inset of Fig. [4|). T he 
ambiguity of the discrimination is quantified in §A 2| 

Applying the homodyne projection |x2 = 0)(x2 = 0| on 
the second mode of Q, we are left with 



l^out) 



s"' \V2a) + |-V^a) ±2e-2"'|0) 
2 v'cosh(2a2) + e-2«^ ± 2e-«' 



(5) 



The fidelity of this state with an ideal even cat state of 
size /3 is 



F = sech(/3^) 



cos 



h(V2. 



±e- 



cosh(2a2) + e-2«' ± 2e- 



(6) 



In practice, a valid output is heralded whenever the ho- 
modyne detector records a state whose a;-quadrature lies 
within a window AQ around x = where AQ shall be 
expressed in shot noise units (SNUs) [40 . We keep track 
of such a post-selection window because no practical ho- 
modyning device has enough resolution to truly project 
onto an exact quadrature |<7o)(^o|- Experimentally, such 
a precise projection would not be desirable either for it 
would lead to very small success probabilities. A com- 
promise is therefore to allow a finite post-selection range. 
A full model of this realistic "wide" homodyning is pre- 
sented in Sec. IA4I 



L.o. 
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■ vacuum 

■ even cat (a = 2) 



FIG. 4: (Color online) Setup for the amplification of two ideal 
Schrodinger cats into a larger even cat. The two inputs are 
mixed on a symmetric beam splitter and one of the outputs is 
projected onto an x-quadrature window of width AQ around 
xo = by an otherwise ideal homodyne detector. Inset: Wave 
functions (x|0) (solid curve) and {x\k-^{2)) (dashed curve) of 
the vacuum and of an even cat state of size a = 2, respectively. 
The two functions are best distinguished at x = where their 
overlap is minimized. 

From now on, we shall only consider odd cat inputs 
to the amplification setup. (The next section deals with 
approximations to odd cat inputs, namely PSSVs.) A 
contour plot of the fidelity between an even cat of size 
P and the amplified state from two odd cats of size a 
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is shown in Fig. [5] for a homodyning window of 1 SNU. 
The amphfication factor that was also witnessed in 
earher schemes [30l [311 [35] is recognizable as the slope 
£ where the fidelity is optimized. The bend of this opti- 
mum crest for a < 1 arises from the vacuum component 
which "survives" the post-selection but vanishes from the 
output state |V^out) foi" larger a. This feature is in a sense 
a manifestation of the discreteness of photon numbers for 
weak coherent states. As can be seen from ([5|, the output 
consists of an even cat minus a vacuum. For a ^ 0, this 
"subtraction" of the vacuum component yields a state 
whose 2-photon component has a relatively higher weight 
than in any even cat of size /3 <C 1. The proportionality in 
amplitude between input and amplified cats thus breaks 
down, and it is instead \j3 ^ 1) — of all even cats — that 
exhibits the best fidelity with the output. 



AQ = 1 SNU 




FIG. 5: (Color online) Contour plot of the fidelity between 
an even cat of size /3 and the amplification obtained from two 
odd cats of size a. The post-selection window is set to 1 SNU. 

Just as in Fig. [2] for the PSSV state, it is informative 
to look at the mapping between input and output pa- 
rameters which optimizes fidelity. This is shown in Fig. 
[6] where the optimal input state amplitude is plotted as 
a function of the amplitude /3 of the target state 
For example, if one wants to produce an even cat state of 
size /3 = 2, then an input odd cat state of size a = 1.4 is 
required. The fidelity of the output with = 2)) will 

be nearly 100% and the success probability of the opera- 
tion about 20%. As discussed above, the fiat plateau for 
P <1 corresponds to the range where the vacuum com- 
ponent that filtered through the post-selection becomes 
predominant. The consequence is that the single-photon 
"cat state", \i<i-{a ^ {))) ~ |1), becomes the only input 
to optimize outputs of target size {3 G [0, 1]. 

To assess the robustness of the amplification scheme 
to the homodyning width AQ, Fig. [7| plots the fidelity 
between the output produced from two odd cats \n-{a)) 
and an even cat \nj^{l3 = \/2a)) given different homodyn- 
ing windows. The fidelity curve for exact homodyning 
AQ ^ is also plotted as a reference. It can be seen 
that the amplification is only vulnerable to AQ for small 




FIG. 6: Top: Maximum fidelity (black) and corresponding 
success probability (gray) of the amplified output with respect 
to an even cat state of size /3. Bottom: Amplitude a of the 
input required to obtain the maximum fidelity of the output 
with an even cat state of size (3. The dotted line marks the \/2 
amplification factor. The homodyne post-selection window is 
1 SNU wide. 



inputs. Beyond a target size of /3 ^ 3.5, homodyning 
widths of up to 8 SNUs hardly have any effect on the 
fidelity. From an experimental point of view, this ro- 
bustness of the homodyning post-selection allows one to 
reach higher success probabilities without compromising 
fidelity. For example, by sending in two odd cats of size 
a = 2.5, one every two homodyne measurements will suc- 
cessfully herald an even cat of size P — a/2 x 2.5 ~ 3.5 
with a fidelity of nearly 100%. 

For another perspective on the dependence of fidelity 
on homodyning width, let's consider the amplification of 
two odd cat states of a fixed size a = 1. The decrease 
in fidelity of the output with an even cat of size /3 = >/2 
is traced in Fig. |8] Also shown are three sample Wigner 
profiles of the output which display a clear degradation 
of the negativity as AQ increases. 

The simulations presented above, as well as all other 
numerical results in this article are arrived at by a generic 
method of simulating linear transformation and projec- 
tive measurements of states consisting of Gaussian su- 
perpositions (see Sec. A3). 



III. 



AMPLIFICATION OF APPROXIMATE 
SMALL CATS 



In this section, we consider the more realistic case 
where PSSVs are amplified, i.e., the inputs to Fig. [4] 
are the outputs of Fig. [l] In Fig. [9) the fidelity profile 
with an even cat of size /3 is plotted with respect to the 
input (pure) squeezing for AQ = 1 SNU and T = 95%. 
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FIG. 7: Solid curves: Fidelity between an even cat of size 
13 = ^/2a and the output of mixing two odd cats of size a 
for five homodyning widths AQ (in SNUs). Dashed curves: 
Corresponding success probabihties. (The curve for AQ 
is zero throughout since the probabihty of picking out the 
exact X = quadrature is vanishingly small.) The shading of 
the curves is lightens with larger homodyning windows. 
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FIG. 8: (Color online) Fidelity of the state amplified from 
two cats of size a = 1 with an even cat of size /3 = 
for a homodyning window AQ G ]0, 15] SNU. The decrease 
in fidelity for wider homodyning windows is understandable 
from the increased overlap between the vacuum and the cat 
state, as illustrated in the inset of Fig. [4] The Wigner profile 
of the output is shown for three sample values of AQ at 1, 5, 
and 10 SNUs, respectively. 



The contour lines of PSSV generation are overlaid to vi- 
sualize the amplification, i.e., the shift of the high-fidelity 
area upwards to larger values of P — cf. Fig. [2j 

In order to assess the performance of the amplification 
scheme, we shall set a target even cat state of ampli- 
tude /3 and retrace the quantum circuit to see what in- 
put squeezing is necessary to achieve the highest fidelity 



with \k:_^{P)) at the output. This is shown in Fig. 10 
along with the dependence of the success probability and 
fidelity on /3, as well as the size a of the odd cat that best 
matches the input PSSV. Let's assume, for instance, that 
we want to produce an even cat of size /3 = 1.5. The re- 
quired squeezing for the PSSV will then be around 2.9 
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FIG. 9: (Color online) Contour plot of the fidelity between 
an even cat \k,-{/3)) and the amplified PSSV obtained from 
an anti-squeezing between and 15 dB. The transmission 
of the subtraction beam splitter is 95%. The contour lines of 
PSSV generation fidelity are overlaid to better visualize the 
shift in target amplitude /3 resulting from the amplification. 
The blank vertical stripe for < 0.8 dB is a region where 
numeric underflow is too frequent to produce reliable data. 
(This is because the normalization factor which enters in the 
fidelity is itself inversly proportional to the success probabil- 
ity. The latter tends to negligible values for small squeezing.) 



dB, corresponding to a fidelity of 96.4% with an odd cat 
state of size a ^ 1.0. The output, however, will have a 
fidelity of 92.6% with |/^+(/3 = 1.5)) and the success prob- 
ability of the amplification will be 20% (notwithstanding 
the success probability of 0.6% required to produce the 
"offline" PSSV). 

Note that the solid curves in the first two plots of Fig. 
[To] only start at f3 ~ 1.2. This is because below that 
threshold, the optimal squeezing lies in a numerically un- 
stable region marked the blank band in Fig. [9j 

A first-hand observation to be drawn from Fig. [lO] 
is the similarity it shows with Fig. [6] for small target 
cat sizes. Both scenarios with ideal cats and PSSV in- 
puts start out with a fidelity around 30-35% which then 
increases to over 90% for targets cats of size /3 ~ 1.5. Be- 
yond this point, however, the performance of PSSV am- 
plification starts to degrade, whereas that of ideal cats 
can be pursued to indefinitely large target sizes while 
retaining unit fidelity. One conclusion is therefore that 
the approximation of cat states from PSSVs can be used 
to generate amplified states of sizes up to /3 = 1.5 with 
practically the same fidelity as if one used ideal odd cat 
states as inputs. On the other hand, both Figs. [6] and 10 
exhibit the plateau in optimal input cat size — or in the 
case of PSSVs, effective cat size — which we discussed in 
the previous section. In that region, P G [0, 1], the ideal 
input cat states or PSSVs cease to have any dependence 
on the target size /3 and discrete states, namely single 
photons, become the optimal input state. 

In addition to simulating our own amplification 
scheme, we have overlaid as finer dash-dotted curves the 
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T = 95% AQ = 1 SNU 




FIG. 10: Top: Dependence of the effective size a of the PSSV 
state input on the target size /3 of the output such that the 
fidehty of the output with |/^+(^)) is maximized. The dotted 
hne marks the ^ = amphfication ratio. The subtraction 
beam spUtter for the PSSV generation is set to 95% and the 
homodyning window to 1 SNU. Middle: Pure squeezing re- 
quired of the input in order to achieve the maximal fidehty 
with an even cat of size /3. Bottom: Maximum fidehty obtain- 
able at the output with an even cat of size /3 (dashed black), 
along with the corresponding success probability of the am- 
plification (dashed gray). Also shown is the maximum fidelity 
of the required input PSSV with an odd cat of size a (solid 
black), and the success probability of the PSSV generation 
(solid gray). (In all three graphs, the finer dash-dotted lines 
refer to the results of Lund et al. in 30 , with, in the bottom 
graph, the black and gray shadings representing fidelity and 
success probability, respectively.) 



results of Lund et al, which we shall refer to as the LJRK 
scheme [30]. Instead of using homodyning, they mix the 
heralding arm with a coherent state on a balanced beam 
splitter such that both emerging modes contain at least 
one photon if the heralded mode is amplified. Contrary 
to the homodyne method, the LJRK projection is unam- 
biguous as it does not yield any residual vacuum compo- 
nent, unlike in ([5|. This explains why the fidelity in the 
LJRK scheme remains quasi-ideal for low target sizes /3. 
Beyond P ~ 1.5, however, both methods are comparable 
in terms of fidelity and amplification factor. The homo- 
dyne method offers nonetheless experimental advantages 
over the photon detection of LJRK in that it does not suf- 
fer as much of quantum inefficiency or electronic noise. In 
particular, the need for coincident detection of photons 



in LJRK suppresses the success probability by the square 
of the quantum inefficiency. (The dash-dotted curve in 
Fig. 10 assumes ideal quantum efficiency.) 

If one assumes ideal quantum efficiencies at the detec- 
tors, the two main nuisances to the performance of the 
amplification scheme are (i) the impurity V of squeez- 
ing, and (ii) the width AQ of the homodyning detection. 
To visualize the robustness of the scheme to these two 
factors, let's choose an optimal input squeezing with a 
target even cat of size, say, /3 = 1.5. From Fig. [Toj this 
corresponds to = —2.9 dB for a fidelity of 92.6% and 
a success probability (assuming offline PSSVs) of 20%. 
From this reference point, the trend of fidelity and suc- 
cess probability varying either V or AQ is plotted in Fig. 
[TT| Recall that, by convention, we model a decrease of 
purity by an increase in anti-squeezing {^x) while main- 
taining squeezing proper (^p) fixed. 

AQ = 1 SNU T = 95% 
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FIG. 11: Robustness of the PSSV amplification to squeezing 
impurity (top) and homodyning width (bottom) . The fidelity 
with an even cat of size /3 = 1.5 and the success probability 
of the amplification (assuming offline PSSVs) are plotted in 
black and gray, respectively. The transmission of the sub- 
traction beam splitter from which the kittens are generated 
is 95% and the anti-squeezing in p is fixed = —2.9. (This 
squeezing is chosen such that the fidelity with \k-\-{/3 = 1.5)) is 
maximized in the case of p ure squeezing and AQ = 1 SNU — 
cf. the middle plot of Fig. [To]) 



IV. SCALABILITY 

As was shown in the previous section, we are bound 
by a trade-off between amplification and fidelity. I.e., 
if the output fidelity is to be kept high, one has to work 
with smaller kittens to begin with, and thus cannot reach 
higher amplitudes (cf. Fig. 10). A high fidelity can 



however be achieved with larger states if we cascade the 
setups into a complete binary tree. We shall discuss the 
scalability of such a recursive amplification in the case 
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of perfect homodyning (AQ 0) and pure squeezing 
(e = 1), both for ideal input cats and for PSSVs inputs. 

Let's describe each iteration stage k as the amphfica- 
tion of two identical wave functions (pk{xi) and (pk{x2) 
into a larger state (pk-\-i{xi) where mode 2 has been post- 
selected by homodyning. This iteration step is given by 



{50:50 beam splitting} 



V2 J ""'K V2 
{projection on \x2 = 0){x2 = 0|} 



X2 



(Pk^i{xi) (X 



^k 



Xi 

71 



(7) 



where the initial wave function (fo is that of the inputs. 
It can be obtained from the 1- variable analog of ( [AIG ). 
In the case of input PSSVs, this involves the highly un- 
balanced mixing two vacua, one of which is squeezed, 
and — ideally — a one-photon projector (f\i)' I.e., 



r 



•</'*i)(a;2)da;2. 



(8) 



Since we only want to investigate how scaling behaves, 
however, we shall simplify the PSSV by a squeezed single 
photon. 



lim (z?o(^) 



(a:|5|l) 



^/2 



1/4^3/2' 



V2s 



TT 



(9) 



where s = 10^/^^ is the factor by which the quadrature 
phase space is re-scaled as a consequence of the pure 
squeezing (i.e., x ^ sx and p ^ p/s). 

The wave function of the state at the k^^ iteration is 
given in closed-form by 



where 



(-) = j^d' (2-/^.) 



^0 



dx 



(10) 



(11) 



is the normalization factor. 
Fig 



12 tracks both the amplitude and the fidelity of 



the outputs as a function of the number of iterations 
k. As already expected from (10), the output amplitude 



grows as >/2 . This does not imply however that itera- 
tions could be carried out indefinitely. For the particular 
choice of initial conditions plotted, a = 1 for the cat or 
3 dB of squeezing for the PSSV, the fidelity drops below 



90% at the fifth iteration. This shortcoming of recursive 
amplifications is due to the non-unitarity of the amplifi- 
cation in the coherent state basis. Looking back at the 
idealized case of (§, we see that the amplification is not a 
straightforward mapping of a cat state onto a larger one, 
but instead introduces an extra vacuum term. This vac- 
uum, which arises from the intrinsic ambiguity of the ho- 
modyne projection, is also amplified along with the cat. 
Even if one does start the first iteration with an ideal 
cat, any subsequent iteration k will inherit this vacuum 
component, which in turn will contaminate the following 
step k -\- 1 with additional terms orthogonal to an even 
cat. 
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FIG. 12: Maximum fidelity (squares) and corresponding ef- 
fective size /3 (disks) of the output as a function of the number 
of iterations. The inputs to the first iteration are an ideal odd 
cat of size a = 1 (black) or a PSSV squeezed by 3 dB (gray). 

The success rate of iterated amplifications is the main 
issue facing scalability. As shown in Fig. [l3j the scheme 
exhibits a sharp drop in the success probability with the 
increasing number k of iterations. Another obvious over- 
head is the number of input states which grows as 2^. 
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FIG. 13: (Color online) Success probability of iterated am- 
plifications for an input cat state \hi-{a = 1)) as well as for 
various input PSSVs with squeezing ^ G {0.5, 1, 2, 3, 6, 9} dB. 
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CONCLUSION 



periods using a particularly reliable locking system. 



We have presented an amplification protocol for cat 
states that is based on imprecise homodyne measure- 
ment. The performance of the scheme was assessed in 
terms of fidelity and success rate and illustrated by an 
optimized relation between the target size f3 and the in- 
put size a of the cat states involved. Given that ideal 
cat states are challenging to produce, we also presented 
how the amplification behaves with approximations to 
cat states, namely photon-subtracted squeezed vacuum. 
Here again, we determined the optimal relation between 
the input squeezing and the effective size of the output. 
The purity of squeezing at the input was determined to 
be crucial in achieving a high fidelity at the output. The 
amplification was however relatively robust to imprecise 
homodyne thresholding, thereby allowing an increase in 
success probability. 

The recursive application of the amplification protocol 
is then simulated in the idealized case of exact homo- 
dyning and pure squeezing. Although the amplification 

/—k 

factor does grow as as v2 with the number k of re- 
cursions, the output fidelity eventually degrades due to 
the non-unitarity of the amplification. One is therefore 
constrained to a finite number of recursions where the 
state is amplified while at the same time retaining a high 
fidelity. 

We saw in ^III that target cat sizes (3 = 1.5 can be 



reached with fidelities up to 93% if one uses PSSVs ob- 
tained from a 5% tapped-off squeezed vacuum. (If one 
assumes ideally squeezed single photons, that fidelity can 
even increase to 98%, cf. the first iteration in Fig. [T2j 
For higher iterations of the amplification, amplitudes of 
/3 ~ 3 can be obtained while maintaining fidelity around 
98%.) Although such approximate states may not per- 
mit fault-tolerant quantum computation, they nonethe- 
less allow for proof-of-principle experiments that require 
effective cat sizes larger than those of basic PSSVs. For 
instance, our amplification protocol can find uses in tele- 
portation [36] or some demonstrations of quantum gates 
O [37] . The question of whether the amplified states can 
be harnessed for any particular use in quantum informa- 
tion will require a feasibility study of its own that pays 
particular attention to the tradeoff between target sizes 
and fidelities. 

Let's conclude with a final note on the practical chal- 
lenges to our protocol. In addition to the issues proper 
to PSSV preparation (and already discussed in [T9l - 
[2T| [24l [25^), the key challenge facing the amplification 
protocol is phase stability. This arises because the pairs 
of interfering PSSVs need to be synchronous, thereby 
leading to very small success rates. This is to be factored 
on top of the already small success rate of an otherwise 
ideal post-selection (Fig. [l3| ). Small success rates are 
compensated for by running the experiment over a larger 
batch of input states. However, this in turns requires 
that the relative phases — of the two interfering PSSVs 
and of the local oscillator — be kept stable for protracted 



Appendix A: Mathematical derivations 

1. Expression for the purity of squeezed vacuum 

Let the variance of the vacuum phase space distribu- 
tion be labeled hy Vq. As will be discussed more for- 
mally in §A3[ the squeezing operation along a quadra- 
ture q consists of a re-scaling of phase space q Sqq 
where Sq = y/Vq/Vo such that the new variance along 
q is Vq. The relationship between the dimensionless re- 
scaling factor Sq and the squeezing ^q (in decibels) is 
given by 

^q = 10 logio 77 = 20 logio Sq^Sq = 10^^/^^. (Al) 

The purity of a state of Wigner function W is given by 
V = 27r jjw^ dxdp, (A2) 

and the Wigner function of a squeezed vacuum state S\{)) 
is 



-(a:/s^)^-(p/sp)^ 



(A3) 



The purity of squeezed vacuum can therefore be shown 
to be 



V 



(A4) 



2. Error in the discrimination between the vacuum 
and a cat state with a homodyne projector 

As mentioned in ^Ilj the purpose of the homodyne mea- 
surement is to collapse the output state onto an even cat 
upon the detection of vacuum; cf. Eq. Q and the inset 
of Fig. [4j The projector for this measurement is 



n 



hd 
|0) 



\x){x\ dx. 



(A5) 



Due to the intrinsic overlap of the wave functions of |0) 
and |A>:+(a)), II|Jj| can only act as an approximate discrim- 
inator between them. In addition, the finite width of the 
quadrature-selection window introduces an approxima- 
tion of its own. The error in the discrimination, namely 
the probability of mistaking an even cat for the vacuum 
is given by 



P(|0))' 



(A6) 
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where, P(|7)) = IV {nhd|7){7|} and I7) G 

{|0), Fig. [m] shows that the discrimina- 

tion is best achieved for smah AQ and large /3. The 
asymptotic convergence to Perr = 0.5 at large AQ or 
small 13 indicates a complete lack of discrimination 
between the states: They become equally likely to be 
inferred by the homodyne projector. 




(3 = |3 = 0.5 ■ 



(3= 1.5 ■ 



-P = 2 



FIG. 14: (Color online) Probability Perr of mistaking an even 
cat state of size /3 for a vacuum state as a function of the 
quadrature acceptance window AQ. 



3. Linear transformations and measurements with 
Gaussian states and operations 

In this section, we shall present the calculational tools 
we have used to perform our simulations. All the quan- 
tum states and operations involved in this article are 
made up of Gaussian superpositions in the Wigner pic- 
ture. I.e., any state p or measurement operator II can be 
written in the Wigner picture as a function 



(A7) 



where G is a Gaussian function in both quadratures 



G{x,p) = go exp {gix'^ + g2X + gsxp + g^p^ + g^p + ^e) 

(A8) 

with constant coefficients g^ G C : k G {0, • • • , 6} and 
^1,^4 G M~. For example, a cat state is made up of four 
such Gaussians 



Wi 



{x,p) = ^G^f2(x,p,7), 



(A9) 



where 



2n (-e-27^ T l) ' 

g-£!;^-p^-2\/2i7p 

~2n (-e-27^Tl)' 
27r T l) 

27r T l) 



whereas the APD operation n^^^^) 
up of only two Gaussians 



(AlOa) 
(AlOb) 
(AlOc) 
(AlOd) 
i- |0)(0| is made 



I^n(APD)(^,p) = '^Gf^^APB,j){x,p) 



where 



Gn(APD,i)(^,p) 
Gn(APD,2) {x,p) 



2^' 



(All) 

(A12a) 
(A12b) 



If we now assume that all transformations underwent 
by the state p prior to measurement are linear, as is 
the case in the setups presented in this article, then the 
generic Gaussian form (A8) is conserved and only the 
quadrature vector 



• • • ,xn,PnY 



(A13) 



of the N modes involved is transformed by a linear map- 
ping 



(A14) 



where the 2A^ x 27V matrix A is determined by the linear 
assemblage of passive elements making up the circuit, 
e.g., beam splitters, phase shifters, squeezers, etc. 

Consider the generic quantum circuit depicted in Fig. 
Tsjwhere once again w e only ha ve Gaussian states and op- 
erators as defined by (A7l|A8 ). The effect of a re-scaled 
quadrature space (A14) on a given Wigner function is 



equivalent to leaving the space unchanged while submit- 
ting the function to the inverse transformation. We shall 
label the transformed multimode Wigner function by a 
tilde such that the mapping (A14) yields 



W{q) W{A-^q) = W{q). 



(A15) 



If we now apply the measurement operators over M < 
N modes, the output over the remaining N—M modes is 
given by 



Pout 



. . AT M 

= yy n ^irJ • n K""^ ^^^dp^ (ai6) 



10 



A 

transform 







CO 
CD 

o 



FIG. 15: (Color online) Sketch of a black box quantum circuit 
made up of A/" = 5 input modes, M — A measurement modes, 
and N—M = 1 output modes. Each state occupying mode k is 
represented by a Wigner function W^^\ The input [In 
is mapped according to ( A15 ) into a linearly transformed state 
Y\n ^\n^ which is then projected onto a "measurement state" 



i Im n 



where the superscript denotes the k^^ mode. 

Note that the number of output modes equals the num- 
ber of input modes N minus the number of measured 
modes M. If M = TV, then there is no output state and 



the inner product ( A16) leaves us with a scalar represent- 
ing the success probability — or norm — of the projection 
of input states onto the measurement operators. Alter- 
natively, this number could be interpreted as the fidelity 
between the overall states at either side of the diagram 
in Fig. [I5l 



The rationale for (A16) is valid regardless of whether 
the Wigner functions are Gaussian superpositions: From 
a mathematical standpoint, projective measurements are 
inner products between the measured state and the mea- 
suring operator. In functional analysis, this translates 
to an integral of the product of two states — the mea- 
sured and the measuring state — over the entire phase 
space of measurement. The simplicity that comes from 



which, in the Wigner representation would most appro- 
priately be given by the Dirac delta function 

WYiT,{q^)=5{q-q^). (A18) 
If we choose for simplicity that = Xq, then the projec- 
tion operation ( A16[ ) onto some state W{x^p) yields 



JJw{x,p)- WuB{qo)dxdp 



p) • 5{x — xq) dx dp 



W{xo,p) 



(A19) 



Intuitively, the probability of measuring exactly the 
quadrature qo should be vanishingly small. Yet, P{xo) 
has a finite value since it is a probability density. To 
model the more realistic scenario where homodyne de- 
tection has a given resolution bandwidth Aq around the 
measured value qo, an interval projector is more appro- 
priate 



nHD(^o,Ag)= / \q){q\dq (A20) 

or, in Wigner representation 



using (A16) with Gaussian superpositions is due to the 






fact that products and integrals of Gaussians are also still 




/ 


Gaussian. 






1 






2^ 





5{q-q')dq' 



Aq]-Q[q-qQ--Aq 



4. Wide homodyne measurement 

A homodyne measurement corresponds to a projection 
on a given quadrature |g'o)(^o| where |g'o) = cos^|xo) + 
sin6>|po)- This projector, when applied to a state whose 
Wigner function is p), "picks out" the cross-section 
along go to yield a probability density. Such a projection 
is given in the quadrature basis by 



(A21) 



where B is the Heaviside step function. Note how (A21) 



nHD(<?o) = ko)(^o| 



(A17) 



does not fit the Gaussian form of ( A8 ) . For calculational 
ease, the wide homodyning operation is therefore best 
performed last so as to maintain the convenient Gaus- 
sianity of the quantum states as long as possible. 

In the limiting case of ideal homodyning, AQ ^ 0, we 
expect (A18) and (A21) to yield the same result. 
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